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ABSTRACT 

We report the discovery of a very low mass binary system (primary mass <0.1 M©) with a projected sepa- 
ration of ^5100 AU, more than twice that of the widest previously known system. A spectrum covering the 
1-2.5 jim. wavelength interval at R ^1700 is presented for each component. Analysis of the spectra indicates 
spectral types of M6.5V and M8V, and the photometric distance of the system is ~62 pc. Given that previous 
studies have established that no more than 1 % of very low mass binary systems have orbits larger than 20 AU, 
the existence of such a wide system has a bearing on very low mass star formation and evolution models. 
Subject headings: binaries: general — stars: low-mass, brown dwarfs — stars: individual (2MASS 
J012655.49-502238.8, 2MASS J012702.83-502321.1) 



1. INTRODUCTION 

A large fraction of the stars in the Galaxy reside in multiple 
systems. Generally, higher multiplicity fractions are found as 
the stellar mass increases, with a binarity fraction reaching 
--80% for B stars (Kouwenhoven et al. 2005). At the other 
end of the distribution, surveys of very low mass (VLM, mass 
< O.IM0) stars and brown dwarfs (BDs), both in the field 
and in star forming regions, have found a binary fraction of 
^20% (Burgasseret al. 2007); corrections for incompleteness 
of tight systems could revise this number to a higher value. 

There are important differences between the VLM binary 
systems and their more massive counterparts, possibly reflect- 
ing differences in their formation mechanism. The mass ratio 
distribution of VLM binaries is strongly peaked at unity (Bur- 
gasser et al. 2007), whereas it is more evenly distributed be- 
tween and 1 for more massive stars. The results of Duquen- 
noy & Mayor (1991) indicate that around 50% of solar-type 
binary systems have separations larger than ^40 AU, and that 
these systems can have separations as large as ^20 000 AU. 
On the other hand, small separations are strongly favored 
among VLM binary systems, less than 1 % of which have sep- 
arations larger than 20 AU (Burgasser et al. 2007). 

To date, only five VLM systems with separation > 
50 AU are known, two in the field, Koenigstuhl lAB 
(1800 AU) and DENIS J055146.0-443412.2AB (220 AU) 
(Billeres et al. 2005; Caballero 2007), and three in star form- 
ing regions: 2MASS J1623361-240221 (212 AU), 2MASS 
J1101192-773238AB (242 AU), and 2MASS J1622252- 
240514 (243 AU) (Jayawardhana & Ivanov 2006; Close et al. 
2006; Luhman 2004). A possible binary system having a sep- 
aration ~1700 AU was recently identified in the a Orionis 
cluster (Caballero et al. 2006), but this system lacks proper 
motion (PM) confirmation of physical association. In this 
Letter we report the identification of 2MASS J012655.49- 
502238.8 and 2MASS J012702.83-502321.1 as part of a 
VLM binary system (hereafter referred to as 2M0126AB) 
with a projected separation of ^5100 AU. 
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2. DISCOVERY OF THE SYSTEM AND OBSERVATIONS 

The system was found by comparing the positions of the 
2MASS point source catalogue (Cutri et al. 2003) with those 
of the /-band Digitized Sky Survey images to search for rel- 
atively high PM (>0."1 yr"') BD candidates at high galactic 
latitudes {\b\ > 30°). Among the sample of candidates found, 
the two sources mentioned above appeared to share the same 
PM. This was further confirmed using /-band observations 
obtained on 2006 June 21 with the Observatoire du Mont Me- 
gantic wide field near-infrared camera CPAPIR (Artigau et 
al., in preparation), currently installed at the CTIO 1 .5-m tele- 
scope. The near-infrared colors (Table 1) of both components 
are indicative of late-M dwarfs (Leggett et al. 2002) and their 
apparent magnitudes place them at a similar distance, suggest- 
ing a physical association (see §4). 

We obtained near-infrared spectroscopic observations 
(R-1700) of 2MASS J012655.49-502238.8 (component A) 
and 2MASS J012702.83-502321.1 (component B) with the 
GNIRS spectrograph (EHas et al. 2006) at the Gemini South 
telescope on 2006 August 11 (see Figure 1) with a 0"3 slit 
under a 0."7-l."0 seeing. For each component, 15 individual 
spectra were obtained for a total integration time of 15 min- 
utes, and their median combination has a signal-to-noise ratio 
of about 30 per resolution element (2 pixels at R^ 1700) in 
the J, H, and K bands. The AO star HD8341 was observed for 
telluric absorption calibration. 

3. PROPERTIES OF THE COMPONENTS 

A summary of all measurements and derived parameters 
is presented in Table 1. The proper motion of each com- 
ponent was calculated from their positions in the SuperCos- 
mos Sky Survey (SSS) /-band catalog (Hambly et al. 2001), 
the 2MASS PSC, and the CPAPIR observations; the coiTe- 
sponding epochs are 1983.661, 1999.778, and 2006.469. A 
correction of (-1-0 ."18,-0 ."10) was appHed to the SSS coor- 
dinates to account for a local systematic offset with respect 
to the 2MASS PSC, as determined from the comparison of 

200 stars located within 15' from 2M0126AB. The astrom- 
etry of the CPAPIR observations was calibrated by using the 
2MASS PSC as a reference frame, and the relative astrometric 
uncertainty after calibration was found to be 0" 16. The uncer- 
tainties on the SSS coordinates (0"15, 0"3) and the 2MASS 
coordinates (0"07,0."07) were estimated from Hambly et al. 
(2001) and the 2MASS PSC, respectively. An eiTor-weighted 
linear regression was performed to obtain the PMs and the 
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TABLE 1 
Parameters of 2MASS J0126AB 
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Fig. 1 . — Bottom panel shows the spectrum of component A with the spec- 
trum of an M6.5 dwarf for comparison. The top panel shows the spectrum 
of B with the spectrum of an M8 dwarf. All spectra are normalized to their 
mean in the range 1.2-1.3 /im and the spectra of comparison dwarfs have 
been offset by +0.4. Spectral regions strongly affected by telluric absorption 
are masked out. The spike-like structure at 2.02 fim in the spectrum of both 
A and B is most likely due to an imperfect correction of a narrow telluric ab- 
sorption band. The comparison spectra are from Rayner et al. (in preparation) 
and were obtained with SpeX at the IRTF. 



corresponding uncertainties. The PM values found for A and 
B, respectively 131 ±9 and 135 ±9 mas/yr in right ascension 
and -53 ± 15 and -47 ±15 mas/yr in declination, agree within 
la. Although the available data make it difficult to quantify, 
one may assume that part of the uncertainty on the individual 
PMs arises from calibration errors common to both objects, so 
that the uncertainties on the difference in PM may be smaller 
than indicated by the above values. 

The spectral types were derived primarily from cross- 
correlation of both spectra with those of known M5-M9 
dwarfs obtained with SpeX at the IRTF'* (Rayner et al., in 
preparation). The wavelength intervals used for the cross- 
correlation were: 1.0- 1.12, 1.18- 1.35, 1.48- 1.8, and 
1 .95-2.4 fj,m. Spectral types of M6.5-M7 and M8 were found 
for components A and B respectively (see Figure 1). The 
spectral indices defined in McLean et al. (2003) yield spec- 
tral types of M6.5-M9 and M6.5-L0.5 for A and B respec- 
tively, consistent with the previous determination. The rela- 
tively short wavelength intervals (0.004 /im) used for calcu- 
lating these indices combined with the modest signal-to-noise 
ratio of our spectra explain the larger uncertainty on the spec- 
tral types using this scheme. The equivalent widths (EWs) of 
the KI lines at 1.169, 1.177, 1.244, and 1.253 ^im for com- 
ponents A and B (see Table 1) are consistent with those of 
M6-M6.5 and M7-M9 dwarfs respectively, based on Figure 
14 and 15 of McLean et al. (2003). Considering all of these 
measurements, we adopt a spectral type of M6.5V for A and 
M8V for B with an uncertainty of ±0.5. 

Photometric distances of 63 ± 5 pc and 61 ± 6 pc are calcu- 
lated for A and B respectively from comparison of their JHKs 
magnitudes with the absolute magnitudes of other dwarfs of 
the same spectral types with known trigonometric parallaxes. 

http://rrtfweb.ifa.hawaii.edu/ spex/spexlibrary/ERTFhbrary.html. 



Parameter 


A 


B 


2MASS designation 


J012655.49-502238.8 


J012702.83-502321.1 


Angular separation 


81f'86±0f'10 


fla COS 


131 ±9 mas/yr 


135 ±9 mas/yr 


Hi 


-53±15 mas/yr 


-47±15 mas/yr 


yb 


21.8 


22.2 




17.2±0.3 


17.6±0.3 


/" 


14.61±0.04 


14.81 ±0.05 




14.05±0.05 


14.16±0.04 




13.68±0.05 


13.62±0.05 


FeH 1.20 nm" 


6.0±1.8 A 


14.5±1.9A 


KI 1.25/im« 


9.8±1.1 A 


10.3 ±1.0 A 


Spectral type 


M6.5V ± 0.5 


M8V ± 0.5 


Photometric distance^ 


63±5 pc 


61±6 pc 


Physical separation 


5100±400 AU 




2670 ± 180 K 


2490 ± 180 K 


Mass (> 1 Gyr)'' 


0.095 ± 0.005 M0 


0.092 ± 0.005 M0 


Mass (~30 Myr)'> 


0.020 ±0.003 Mq 


0.019 ±0.003 Mq 



"Calculated using the SuperCosmos Sky Survey /-band positions (epoch 
1983.661), the 2MASS point source catalogue position (epoch 1999.778), 
and our measurements made with CPAPIR (epoch 2006.469). 

''Assuming V -/ = 4.6 for both objects (see Dahn et al. 2002, fig 3) 

'^From the SuperCosmos Sky Survey catalogue. 

"■From the 2MASS point source catalogue. 

^Spectral indices defined in Gorlova et al. (2003) 

^Determined by comparing the /, H, and Ks magnitudes of A and B with 
the absolute magnitudes of eight M6.5V-M7V and four M8V dwarfs with 
known trigonometric parallaxes; the quoted uncertainties correspond to the 
dispersions of the distances obtained. 

^Spectral-type - 7^// relation of GoUmowski et al. (2004). 

''Determined by minimizing the sum of the square difference between the 
absolute J, H, and Ks magnitudes of 2M0126AB and evolution models 
(Chabrier et al. 2000). The quoted uncertainties are derived from uncertain- 
ties in the photometric distance. 



Using the spectral type-effective temperature (Teff) relation 
of Golimowski et al. (2004), we infer T,ff = 2670 ± 180 K 
for A and T^/f = 2490 ± 180 K for B. According to evolu- 
tion models (Chabrier et al. 2000), such temperatures imply 
masses below 0.1 Mq, irrespective of the system age. 

4. PROBABn.ITY OF RANDOM ALIGNMENT 

Assuming a late-M dwarf spatial density of 2.2 x 10"^ pc~^ 
per /-band magnitude interval (Phan-Bao et al. 2003), there 
should be ~2700M6-M8 dwarfs (12.5 <M/ < 14.5) between 
55 and 68 pc from the Sun. Considering that an 82" ra- 
dius disk occupies 3.95 x 10"** of the sky, and assuming an 
isotropic distribution of late-Ms out to 68 pc, the probability 
of finding at least one such pair of objects over the whole sky 
from chance alignment is 0.15, regardless of their apparent 
proper motion. A Monte Carlo calculation was performed to 
determine the probability that two M dwarfs located in the 
direction of 2M0126AB have PMs larger than 0."1 yr"' that 
agree with each other within our uncertainties at the 2 a level. 
Using the observed dispersions in (U, V, W) of late-M dwarfs 
in the Solar neighborhood (Bochanski et al. 2005), including 
the disk and halo components, and converting these veloci- 
ties into PMs at the position of 2M0126AB, the probability is 
1 .4 X 10"^. Thus, the combined probabiUty of finding a pair of 
late-M dwarfs within 82" of each other and sharing the same 
proper motion within our uncertainties is 2x10"^. Therefore, 
2M0126AB is very likely a bound system rather than a ran- 
dom alignment of two unrelated objects. 
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5. DISCUSSION 

As shown in Gorlova et al. (2003), the EWs of FeH and 
KI provide a rough constraint on surface gravity. Indeed, the 
strength of these features in late-type M dwarfs is systemat- 
ically higher (by a factor 1.5-2.0) for field objects compared 
to young (< 10 Myr) cluster members. Based on Figure 8 
of Gorlova et al. (2003), the measured FeH 1.20/im and KI 
1.25/im EWs of components A and B are in better agree- 
ment with those of relatively old field objects. The mass of 
each component can be determined from comparison of near- 
infrared absolute magnitodes with evolution models (Chabrier 
et al. 2000) assuming an age for the system. If we suppose 
an age greater than 1 Gyr, as suggested by the gravity sensi- 
tive spectral indices, we derive masses of 0.095 ± 0.005 and 
0.092 ± 0.005 M0 for the two components. 

Although the FeH and KI EWs are consistent with rela- 
tively old field objects, the age of 2M0126AB remains poorly 
constrained given the small difference in EW between field 
and very young dwarfs and the significant dispersion of those 
values (see Figure 8 of Gorlova et al. 2003). In fact, there 
is some evidence that 2M0126AB might be part of a young 
association. On the sky, it lies in the core of the young 
(~30 Myr) Tucana/Horologium (TH) association (Zuckerman 
& Song 2004 and references therein) located at a distance of 
37-66 pc, consistent with that estimated for 2M0126AB. Fur- 
thermore, both the amplitude and direction of its PM are sim- 
ilar to those of nearby members (see Figure 2). We note, how- 
ever, that the amphtude of the PM of 2M0 1 26AB ( 1 42 mas/yr) 
is slightly larger than the largest PM of nearby members 
(124 mas/yr). The PM dispersion of members of the asso- 
ciation is ^ 1 1 mas/yr, and is representative of relatively mas- 
sive members (>0.5 Mq). The dispersion is expected to be 
significantly higher for lower mass members (see Figure 4 of 
Pinfield, Jameson & Hodgkin 1998). An age of 30 Myr for 
2M0126AB would imply that the two components are BDs 
with masses of 0.020 ±0.003 and 0.019 ±0.003 M©. Future 
determinations of the trigonometric parallax, the radial veloc- 
ity, the hthium abundance, the surface gravity, and the Ha or 
X-ray emission of 2M0126AB will help in settling this issue, 
and in refining the physical characterization of the system. 

In the scenario where 2M0126AB is a member of the 
TH association, the random alignment probability determined 
earlier needs to be re-evaluated; the major argument in favor 
of a physical association now being the relatively small angu- 
lar separation of the system given the low stellar density of the 
association. The association has a diameter of approximately 
30 pc, hence a volume of 1.4 x 10"* pc^. The uncertainties on 
the distance of both components and their angular separation 
place them within a volume 5.7 x 10^ times smaller. Assum- 
ing a Kroupa three-segment IMF (Kroupa 2001), and know- 
ing that there are 42 members of the TH association between 
B8 and K5 (roughly 4 Mq to 0.67 M©), we estimate that the 
association comprises 32 objects in the mass interval between 
0.015 Mq and 0.025 Mq. Given these numbers, the probabil- 
ity that the two components are not physically bound is only 
8.7 X 10""*. Even assuming a factor of 5 underestimate of the 
number of low-mass objects in the association, this probabil- 
ity would still be only 2 %. 

At a distance of 62 pc, the measured angular separation of 
the system corresponds to a projected physical separation of 
~5 100 AU, more than twice that of the Koenigstuhl 1 AB sys- 
tem (Caballero 2007) and ^20 times greater than those of 
other known wide VLM binaries (see Figure 3). Because 
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Fig. 2. — (a) The PM of all TH member.s listed in Zuckerman & Song 
(2004) and of 2M0126AB are indicated by an arrow originating from their 
right ascension and declination. The position of 2M0126AB is indicated by 
a cross. 2M0I26AB lies in the core of TH and its tangential motion is very 
similar to that of nearby members, (b) Position and PM of 2M0126AB and of 
all stars in the Hipparcos catalogue located within 15 degrees of 2M0126AB 
and with 100 < PM (mas/yr) < 150. The position of 2M0126AB is indicated 
by a cross. There is no overall tendency for stars in this part of the sky to 
move in a direction similar to that of 2M0126AB, or of members of TH. 



of projection effects, the real separation is expected to be, 
on average, 1.4 times larger (Couteau 1960) which, in the 
present case, would be ~7000 AU. Prior to the discovery of 
2M0126AB, projected separations of ~5000 AU had been ob- 
served only for systems of total mass > 1 Mq. The low total 
mass and large separation of 2M0126AB raise the question 
of its dynamical stability. An order of magnitude estimate 
of the expected survival time of 2M0126AB can be obtained 
from the work of Weinberg et al. (1987). Stellar and giant 
molecular cloud (GMC) encounters are the dominant source 
of disruption for low-mass binaries, and these contributions 
are dependent on the mass/separation ratio only, implying that 
the results derived for more massive binaries can be scaled 
to the 2M0126AB system. Given this scaling relation and 
Figure 2 of Weinberg et al. (1987), the stellar and GMC en- 
counter disruption timescale of 2M0126AB is ~200 Myr if it 
is a member of the 30 Myr old TH association or '--^1 Gyr if 
it is a field object. The first timescale above is well beyond 
the age of the TH association, while the second may indicate 
that 2M0126AB is likely younger than 2-3 Gyr. Further, the 
second timescale above is long enough for a significant num- 
ber of similar binaries produced throughout the history of the 
Galaxy to have survived to this day. 

The formation process of VLM stars and BDs is the sub- 
ject of an active debate, partly because their mass is less than 
~0.7 Mq, the typical Jeans mass in molecular clouds (Meyer 
et al. 2000). However, it has been shown recently (Boyd & 
Whitworth 2005; Padoan & Nordlund 2004; Padoan et al. 
2005) that turbulence can lead to fragmentation of unstable 
cores with masses as small as ~0.003 Mq. A possibihty is 
that VLM stars and BDs are such low-mass embryos dynam- 
ically ejected out of their feeding zone during the early stages 
of formation (Reipurth & Clarke 2001). This is a popular sce- 
nario because it agrees with the observation that most VLM 
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Fig. 3. — Orbital separation vs total mass of binary systems. VLM systems 
are shown with empty circles and stellar systems are shown with crosses. The 
left and right filled circles show 2M0126AB if it is aged 30 Myr (i.e. a mem- 
ber of TH) or > 1 Gyr respectively. The labels refer to the large separation 
VLM systems mentioned in the text. 



binaries have close separations, weakly bound wider binaries 
being easily disrupted by the ejection. Nevertheless, detailed 
simulations (Bate & Bonnell 2005) of this mechanism have 
succeeded in producing some VLM binaries with separations 
~100 AU, similar to the four previously known VLM bina- 
ries mentioned above albeit with lower mass ratios. These 
wide binaries result when low mass cores become bound fol- 
lowing ejection from the cluster at the same time and in the 
same direction. It remains to be seen whether binaries Uke 
2M0126AB, with orbital separations ~20 times larger, can be 
produced by this mechanism. 

Photo-evaporation of pre-stellar cores by the strong ioniz- 
ing radiation of nearby massive stars (Whitworth & Zinnecker 
2004) is another possible VLM star and BD formation pro- 
cess. However, the presence of massive stars required for this 
to happen would most likely disrupt a VLM binary such as 
the one described here. Formation of the secondary by direct 
gravitational collapse within the accretion disk of the primary 



(Boss 2000, 2006) is yet another possibility, but it is unlikely 
to form nearly equal mass binaries and would require a disk 
of unreasonably large spatial extent. 

The formation of BDs following the encounter of proto- 
stellar disks has been shown possible in simulations (Shen 
& Wadsley 2006). In this case, the disks are normally sta- 
ble against gravitational collapse; it is only the violent col- 
Usion between two disks that triggers instabilities leading to 
the formation of BDs. Although current simulations are not 
sufficient to address directly the question of wide binary BD 
formation, the results indicate that about one third of the BDs 
produced are ejected from the system, and wide binaries may 
result as in the ejection scenario discussed above. So far, 
the simulations have only produced BD mass objects and not 
VLM stars through this mechanism. 
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